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Abstract A;— point correlations of complex zeros for Gaussian ensembles of Random Poly- 
nomials of order N with Real Coefficients (GRPRC) are calculated exactly, following 
an approach of Hannay |J for the case of Gaussian Random Polynomials with Com- 
plex Coefficients (GRPCC). It is shown that in the thermodynamic limit N —>■ oo of 
Gaussian random holomorphic functions all the statistics converge to the their GRPCC 
counterparts as one moves off the real axis, while close to the real axis the two cases are 
essentially different. Special emphasis is given to 1 and 2 point correlation functions in 
various regimes. 



The problem of statistics of zeros of random polynomials of order N, and of random 
holomorphic functions as N — > oo in general, arises in various contexts in quantum chaos 
|| |||. The motivation for this work was the problem of statistics of zeros of coherent state 
(Husimi) or Bargmann Q representation of eigenstates of chaotic systems ^ It has 
been conjectured || that zeros of Bargmann or Husimi representation of an eigenfunction 
of 1-dim classically chaotic system should be uniformly and randomly scattered over the 
classically chaotic region of phase space. Bargmann representation of an eigenstate is an 
entire analytic function in a complex phase space variable z = q + ip, sometimes it is 
even a polynomial of a finite order, like for example in the case of spin systems where 
the phase space manifold is a sphere parametrized by (#,</>) and z = cot (0/2) exp(i0) is 
a stereographic projection. The coefficients of a power series of such entire functions or 
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polynomials are just the coefficients of an expansion of the chaotic eigenstate in a com- 
plete set of (say harmonic) wavef unctions. Applying the random matrix theory one argues 
that these coefficients should be uncorrelated (real/complex in the presence/absence of 
anti-unitary symmetry) pseudorandom Gaussian variables. Thus one can introduce the 
statistical ensembles of random polynomials of order iV (or random analytic functions in 
the limit iV — ► oo) and argue that statistical properties of their zeros can be used as a 
model to describe statistical properties of zeros of a Bargman representation of chaotic 
eigenstates of real systems. 

Recently, Hannay [H has calculated general fc-point correlation functions of zeros of a 
random spin state in a coherent state representation which is described by the random 
polynomial with uncorrelated complex Gaussian coefficients, and solved the problem of 
statistics of zeros of GRPCC — Gaussian random polynomials with complex coefficients in 
general. It has been demonstrated numerically |]7|, [| that his results on GRPCC provide 
a universal description of the statistics of zeros of Bargmann or Husimi representation 
of chaotic eigenstates for systems without an anti-unitary symmetry. Here we adopt this 
approach and solve the general problem of statistics of zeros Zk of GRPRC — random 
polynomials f(z) of order N 

TV N 

f(z) = Y, a nZ n = a N Y[(z- Zj ) (1) 

n=0 j=l 

with real (Gaussian) coefficients a n . We argue that the obtained results may be used to 
describe statistics of zeros of eigenstates of 1-dim. (and quantum Poincare sections || 
and other reductions [l^, |ll[] of 2-dim.) chaotic systems in Bargmann representation with 
time reversal invariance (or any other anti-unitary symmetry Q) to the same extent as 
Gaussian orthogonal ensembles of random matrices can be used to describe the Hamilto- 
nian and the typical observables. 

In the literature one may find several results on the distribution of complex zeros of 



random polynomials with either complex |J or real |12| Gaussian coefficients (see also || 
and references therein). The formula for 1-point function given below (|IP|) (in the special 
case where the variances of all coefficients are equal) is equivalent to the theorem 1.1 of 
Shepp and Vanderbei [ |12| . 

Take a A;— tuple of complex numbers z = (z\, . . . , Zk)- Since a n are real Gaussian random 
variables (which in general need not be uncorrelated!), their real linear combinations 



fJ = Ref( Zj ) } f; = Imf( Zj ), f>r = Re-f( Zj ), ft = 1m—f( Zj ), j = 1, 



(2) 



are also real Gaussian random variables with a joint distribution 

P(f , f, f r , f) = (det 27TM)- 1 / 2 exp (-\{f , f , f' r , f) ■ M\f r , f, f' r , f)) . (3) 

2 For a general anti- unitary symmetry, the coefficients of the random polynomials (Q) are of the form 
o-n = r n e ie " where r n are real Gaussian random variables and 9 n are fixed (nonrandom) phases (which 



determine the symmetry curve in complex z-plane, such as in fig. 6 of ref. H). Then one may use the 
same general approach described below, eqs. (&U3). 
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M is a 4k x 4k real symmetric positive correlation matrix 
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where () denotes the Gaussian ensemble averages which can be calculated using (|l],||) in 
terms of input data {a n a m ). One can write the k— point correlation function pk{z) in the 
following form 



(5) 



Pk (z) = / p(o,o,r,r) mfj r ) 2 + (sifWjtfi 

where the factors (/j r ) 2 + (/j 2 ) 2 are just the Jacobians of transformations from the pairs 
of real variables (/J, /j) to complex variables — zeros Zj. The integral can be written in 
terms of derivatives of a generating function Z k {u, v) 



(6) 



which is an ordinary Gaussian integral and can be explicitly calculated 



Z k (u, v) = (det 27TM)- 1 / 2 / exp {-\{f, f) ■ L(f r , f) + if ■ « + if 1 ■ «) II d f'i d f 

J 3=1 

= (det2vrA) _1/2 exp (-\{u,v) ■ L(u,v 



(7) 

where L = C — B A Bisa lower right block of M and we have used an identity || 
det L/ det M = 1/ det A. At this point it is convenient to switch on the equivalent complex 
variables f = f + if 1 , f = f' r + if' 1 , w = u + iv and their complex conjugates. Then 
one can write eq. as 



Pk{z) 



u/2 II dfofttf ex P {-\i w *i w ) ' L (™' w *i) \w=o 



(det2vrA) 1 /2 ±± 



(det2 7 rA)- 1 / 2 n d Wj d w * {{w*,w) ■ L(w,w*)) k \ w=0 



3=1 



where all the 2k x 2k real matrices should be transformed by the rule 

1 1 



X = U f XU, U = i 



i\ -il 



jiving L = C - B'A B with 

A 

B 

C 



(fjfi) 


(fjfk) \ 


\jj Jk j 


(&»)) 


(fsft) 


ML) > 


1 f* fl*\ 
\Jj Jk 1 




(m) 


(f'sfi.) 


(frfk) 


(f?fi) 



At, 
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Applying a little combinatorics on eq. we finally obtain the general result 



sper (C - B* A^B) 
Pui») = V i ^ (12) 



Vdet 2?r A 

where we introduce the semi-permanent of a 2/c x 2/c matrix 



s 



perL = E E II L jr+fc,Jp (r ) ( 13 ) 

Ji< --<7fc p€5 fc r=l 
fi<-<ifc 

The first sum runs over (2k)\/(k\) 2 ordered combinations of k out of 2k indices j m and 
their complements l n while the second sum runs over k\ permutations p of the symmetric 
group Sfc. The sum of indices j r + k should be taken modulo 2k. 

So far we have not assumed anything about the correlations between the coefficients 
a n expect the Gaussian nature of the joint distribution of coefficients a n . Now we shall 
asume that Gaussian coefficients a n are uncorrelated and define the polynomial g(s) with 
positive coefficients b n — the variances of a n 

(a n a m ) = b n 8 nm , b n > 0, (14) 

N 

g(s) = (15) 

n=0 

The matrices A,B, and C can be easily expressed solely in terms of a polynomial g and 
its derivatives g', g" 

M*) = 9(zjz!), (16) 
B jt (z) = d z: A jl (z) = z J g'(z,z*), (17) 
C fl (z) = d Zj d zt A jl (z)=g'(z j zt)+z j ztg"(z j z:) (18) 

where we let indices j, I to run from 1 through 2k and put Zk+j '■— z*. Note that the 
time-reversal symmetry — the symmetry of zeros with respect to the reflection over the 
real axis is present also in the k— point correlation functions, namely 

p k (zi, ...,Zj,..., z k ) = Pk(zi, ...,z*,..., z k ). 

Without loss of generality one may assume that all points Zj lie on the upper complex 
halfplane, Imzj > 0. Otherwise one gets long range correlations in cases where one of the 
points Zj comes close to the mirror image of one of the other points z\. 

In general, only the 1-point function p\{z) — - the density of zeros is simple enough to 
be written out 

M go + \A 2 9o , (z 2 9-9+ + z* 2 g + g'_)g' Q - \z\ 2 (g' + g'_ + g' 2 )g 

Pl[) Agl-g+g-) 1 ' 2 Agl~g + g-f /2 ' 1 j 

where = ^(kl 2 ); g+ = g(z 2 )i g- = g( z * 2 )- Writing z = x + iy and carefully expanding 
for small y one finds 

Pl (z) = h{x 2 )\y\ + 0{y*), y^O (20) 
h(s) = (2n)- l (gg'-sg' 2 + sgg"y 3/2 (2goi2 + 2(2g 013 -g 112 -g 022 )s 

+ (3^122 - 40X13 + g0u)s 2 + (5-024 ~ 5114 ~ 5033 + 2^123 - 5222)s 3 ) 
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where g = g(s), g nm i = g^ n \s)g^ m \s)g^(s). So quite generally, the density of zeros 
decreases linearly as we approach the real axis. To evaluate the density of zeros on a real 
axis y = one should use a different approach described in ||. In another asymptotical 
regime \z\ — ► oo, only the highest power terms of g contribute, and one finds 

wW = *£^2£(i + (i)V (2i, 



So, the density of zeros vanishes asymptotically since the total number of zeros N is finite. 

Now we shall study the thermodynamic limit N — > oo. It is convenient to study ran- 
dom holomorphic functions which provide a uniform distribution of zeros in the complex 
plane. A unique choice (up to rescaling s — > As) is b n = l/n\ giving 

</(*)= exp(s). (22) 

Such random holomorphic functions naturally arise when one studies Bargmann repre- 
sentation of 1-dim chaotic eigenstates in the usual (p, q) G 3? 2 phase space. We argue 
that any other choice will only affect the density of zeros pi(z) while properly rescaled 
local statistics should be independent on the choice of g{s) provided that variances of 
coefficients b n depend smoothly on n. 

Away enough from the real axis InxZj ^> 1 one may neglect the offdiagonal k x k blocks 
of matrices A, B, C since the ratios of the corresponding matrix elements become expo- 
nentially small | exp(zjZi)/ exp(zjZi)\ = exp(— 2lmzjlmzi). Then using straightforward 
results 

2- fc s P er( Ll1 ° T ) = perLn := £ f[ L jAj) (23) 

det( A n A ° T J = (detA n ) 2 , (24) 

where Qn denotes the upper-left k x k block of a 2k x 2k matrix, one arrives to the result 
which is equivalent to the statistics of zeros of GRPCC 



GRPCC/ n per(Cn - B| 1 A 11 1 B 



p k {z)^p^^{z) = ^^^n, ; ag j (25) 

det 7TAn 



To conclude we give some explicit results about 1 and 2 point functions. The density 
of zeros which is shown in figure 1 reads 

/ , • x l-(% 2 + l)exp(-V) 
Pl(X + ^ )= vr(l - exp(-V))3/ 2 ' (26) 

which is a constant 1/tt provided that we are away enough from the real axis. The excess 
of zeros due to the presence of real axis /^(l/V — px(x + iy))dy = 1/n is on the other 
hand just the linear density of real zeros on the real axis! 

The 2-point correlation function P2{zi,z 2 ) is already too lengthy to be written out in 
general. The behaviour of a normalized 2-point correlation function p2{zi, z-i) / p\{z\)j P\{z-i) 
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as we approach the real axis, is shown in figure 2, while far away Imzi } Imz 2 3> 1 it be- 
comes isotropic and the result for GRPCC || applies 



P2(zi,z 2 ) -> <f(\zi - z 2 \ 2 ), 

exp(-2s)(exp(s) - 1 - s) 2 + exp(-s)(exp(-s) - 1 + s) 2 
-k z {\ — exp(— s)) A 

In the asymptotic regime Imzj ^> 1 one can calculate also the number variance Yi 2 {r): 
the variance of the number of zeros J\f(r) inside a circle of radius r 

E 2 (r) = (A/- 2 (r))-(Af(r)> 2 . (28) 

It can be expressed in terms of a four- fold integral (over zi,z 2 ) of a 2-point correlation, 
which can be reduced using eq. (|27|) to a single integral 

£ 2 (r) = r 2 (l — r 2 ) + 87rr 4 J (arccos Vi- y/t(l -t))ip{±r 2 t)dt. (29) 

The number variance S 2 (r) starts as "Poissonian" (Af(r)) = r 2 for small r whereas for 
larger r it has a linear asymptotics (see figure 3) 



oo 
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S 2 (r) = or + 0(l/r) « <TyJ(Af{r)), a =- J s 2 (l - ir 2 (p(s 2 ))ds w 0.36847. (30) 

^ o 

Note that this formula (E9LI30]) is valid also for GRPCC in general. 



In the present paper the statistics of zeros of Gaussian random polynomials with real 
coefficients have been solved analytically (|1^) following an approach of Hannay for the 
case of complex coefficients. Several important special cases have been considered in de- 
tail: (i) the case of mutually uncorrelated coefficients, which corresponds to the Bargmann 
representation of chaotic eigenstates in random matrix regime, has been studied and it has 
been shown that all A;— point correlation functions converge to those of random polyno- 
mials with complex coefficients derived by Hannay as all points Zj,j = 1 . . . k move away 
from the real axis Imzj ^> 1 (^3), (ii) 1-point function - the density of zeros have been 
written out in general (eqs. (p!9| , p6|) and fig. 1) and linear decrease of density towards the 
symmetry line - real axis has been found (|2"0"D, (ii) 2-point function close to the real axis 
have been explored numerically (fig. 2) while simple analytic formula (p7[), which holds 
far away from the real axis (and holds generally in the case of complex coefficients), have 
been used to derive a simple expression for the number variance of zeros inside a circle of 
a given radius (eqs. ( p9| , |30D and fig. 3). 

Discussions with P.Leboeuf, J. H. Hannay, M.Saraceno and K.Zyczkowski as well as the 
hospitality of Institut Henri Poincare (Paris) are gratefully acknowledged. This work 
has been financially supported by the C.I.E.S (France) and the Ministry of Science and 
Technology of the Republic of Slovenia. 
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Figure captions 



Fig.l We show the density of zeros p\(x + iy) in the thermodynamic limit N — > oo given 
by eq. (|26| ) as a function of the distance from the real axis. 

Fig. 2 The normalized 2-point correlation function p 2 (xi+iy, x 2 +iy) / p\{x\+iy) / pi(x 2 +iy) 
in the limit N — > oo between two points, x% + iy and x 2 + iy, which have the same 
distance from the real axis y is shown as a function of \x 2 — xi| for different values of 
y = 0.1,0.3,0.5,0.7,0.9,1.1,1.3,1.5. Note that all curves go to zero as oc y 2 and that 
for y > 1.5 the 2 point correlation function has practically converged to the isotropic 
asymptotic one. 

Fig. 3 The number variance £2(7") i n the asymptotical regime iV — > oo,Imz ^> 1 is shown 
as a function of radius r (|29f) . 
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